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Abstract — In this paper we study rebalancing strategies for 
a mobility-on-demand urban transportation system blending 
customer-driven vehicles with a taxi service. In our system, a 
customer arrives at one of many designated stations and is 
transported to any other designated station, either by driving 
themselves, or by being driven by an employed driver. The 
system allows for one-way trips, so that customers do not have 
to return to their origin. When some origins and destinations 
are more popular than others, vehicles will become unbal- 
anced, accumulating at some stations and becoming depleted 
at others. This problem is addressed by employing rebalancing 
drivers to drive vehicles from the popular destinations to 
the unpopular destinations. However, with this approach the 
rebalancing drivers themselves become unbalanced, and we 
need to "rebalance the rebalancers" by letting them travel back 
to the popular destinations with a customer. Accordingly, in this 
paper we study how to optimally route the rebalancing vehicles 
and drivers so that stability (in terms of boundedness of the 
number of waiting customers) is ensured while minimizing the 
number of rebalancing vehicles traveling in the network and the 
number of rebalancing drivers needed; surprisingly, these two 
objectives are aligned, and one can find the optimal rebalancing 
strategy by solving two decoupled linear programs. Leveraging 
our analysis, we determine the minimum number of drivers 
and minimum number of vehicles needed to ensure stability 
in the system. Interestingly, our simulations suggest that, in 
Euclidean network topologies, one would need between 1/3 and 
1/4 as many drivers as vehicles. 

I. Introduction 

In this paper we study vehicle routing algorithms for a 
novel model of urban transportation system, which involves 
blending customer-driven vehicles with a taxi service. Our 
proposed car- share system is an example of a Mobility- 
on-Demand (MOD) system, and aims at providing urban 
dwellers with the tailored service of a private automobile, 
while utilizing limited urban land more efficiently (e.g., 
by minimizing the automobiles that sit unused) [1]. In our 
system, a customer arrives at one of many designated stations 
and is transported to any other designated station, either 
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by driving themselves, or by being driven by an employed 
driver. The system allows for one way trips, so that customers 
do not have to return to the same stations from which they 
picked up their vehicles. In a typical one way car-share 
system (e.g. Car2Go) it has been observed empirically [2], 
and shown analytically [3], that vehicles become unbalanced, 
accumulating at popular destinations and becoming depleted 
at less popular ones. Our proposed system addresses this 
problem by employing rebalancing drivers to drive vehicles 
from the popular destinations to the unpopular destinations. 
However, with this approach the rebalancing drivers them- 
selves become unbalanced, and hence we need to "rebalance 
the rebalancers" by letting them travel back to the popular 
destinations with a customer. In such a trip, the rebalancing 
driver operates the vehicle as a taxi, driving the customer to 
their desired destination. The system is illustrated in Fig. [T] 
The main difficulty in such a system, and the focus of this 
paper, is how to determine the rebalancing trips and the taxi 
trips in order to minimize wasted trips, while providing the 
best possible customer experience. 

Specifically, the contribution of this paper is twofold: we 
study routing algorithms for the MOD system illustrated 
in Fig. [T] that (1) minimize the number of rebalancing 
vehicles traveling in the network, (2) minimize the number 
of drivers needed, and (3) ensure that the number of waiting 
customers remains bounded. Second, leveraging our analysis, 
we determine the relation between the minimum number 
of drivers needed and the minimum number of vehicles 
needed to ensure stability in the system; these relations would 
provide a system designer with essential structural insights 
to develop business models. Interestingly, our simulations 
suggest that, in Euclidean network topologies, one would 
need between 1/3 and 1/4 as many drivers as vehicles, and 
that this fraction decreases to about 1/5 if one allows up to 
3-4 drivers to take a trip with a customer. 

This paper builds upon the previous work of the authors 
in designing optimal rebalancing policies for MOD systems 
leveraging autonomous operation of the vehicles [4], [3], 
i.e., without the need of human drivers. On the contrary, the 
system proposed in this paper would use technology that is 
available today (i.e., by employing human drivers instead 
of autonomous cars), and our finding are readily applicable 
to existing one-way car-share systems, which already employ 
drivers to rebalance cars using heuristic methods [2]. Further- 
more, by comparing the results in this paper with those in [4], 
one can quantitatively assess the relative benefits of "hi-tech" 
autonomous MOD systems versus "low-tech" driver-based 
MOD systems. The problem addressed in this paper has 



also many characteristics in common with the well-known 
Dynamic Traffic Assignment (DTA) problem [5], [6], [7], 
[8]. The key difference between rebalancing in MOD systems 
and the DTA problem is that in the former the optimization 
is over the empty vehicle trips (i.e., the rebalancing trips) 
rather than the passenger carrying trips. 

The rest of the paper is structured as follows. In Section |TX| 
we present a model for our system with customers, vehicles, 
and drivers represented as a continuous fluid, and we for- 
mally state the problem of rebalancing the vehicles and the 
drivers. In Section [ill] we (i) study the well-posedness of the 
model and characterize its set of equilibria; (ii) determine 
the minimum number of vehicles and drivers needed to meet 
the customer demand; and (iii) show that with rebalancing 
vehicles and drivers the system is indeed locally stable (i.e., 
stable within a neighborhood of the nominal conditions). In 



Section [IV] we show how to optimally route the rebalancing 
vehicles and drivers so that stability (in terms of boundedness 
of the number of waiting customers) is ensured while min- 
imizing the number of rebalancing vehicles traveling in the 
network and the number of rebalancing drivers needed; re- 
markably, these two objectives are aligned, and one can find 
the optimal rebalancing strategy by solving two decoupled 
linear programs. In Section [V] we study the relation between 
the minimum number of drivers needed and the minimum 
number of vehicles needed. Finally, in Section [Vl] we give 



conclusions and discuss future research directions. 




Fig. 1. At each station there are three queues: customers (yellow dots), 
drivers (red dots), and vehicles (small car icons). There are three modes of 
use for a car: A customer can drive a car between stations; a customer can 
be driven between stations by a driver; or, a driver can drive a car between 
stations to rebalance. 



II. Modeling the Mobility-on-Demand System 

In our prior work [3] we proposed a fluid model for 
mobility-on-demand systems and formulated a policy to op- 
timally rebalance vehicles assuming that they could operate 
autonomously. In this paper we consider rebalancing the 
vehicles through the use of dedicated personnel that are 
employed to drive the vehicles. In this section we extend 
the fluid model in [3] to capture the later scenario. 
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TABLE I 

Description of notation for station i 
Definition 



number of customers at station i 

number of vehicles at station i 

number of drivers at station i 

rate of arrival of customers at station i 

departure rate from station i 

travel time from station i to station j 

fraction of customers at station i destined for station j 

rate of rebalancing vehicles from station i to station j 

J2j an 

rate of rebalancing drivers from station i to station j 
fraction of customers traveling from i to j willing 
to use taxis 
Heaviside function 



Basic model: The model in [3] can be formalized as 
follows. Consider a set of n stations, J\f = {l,...,n}, 
defined over an extended geographical area (see Figure [I]). 
Since the model is a fluid approximation, the number of 
customers, vehicles, and drivers are represented by real 
numbers. Customers arrive at station i at a constant rate 
A^ G M >0 . The number of customers at station i at time 
t is Ci(t) G M>o, and the number of vehicles waiting idle 
at station i at time t is Vi(t) G M>o- The total number 
of vehicles in the system is V G M>o- The fraction of 
customers at station i whose destination is station j is p^ 
(where p^ G M>o, Pa = 0, and YjPij = 1). The travel 



time from station i to station j is XL G 



When there 



are both customers and vehicles at station i (i.e., Ci(i) > 
and Vi (t) > 0), then the rate at which customers (and hence 
vehicles) leave station i is when, instead, Ci(t) = but 
Vi(t) > the departure rate is A^. A necessary condition for 
the total number of customers at station i to remain bounded 
is that fii > Xi ; we will assume Hi > Xi throughout the paper 
(the case /^ = A^ can be addressed with techniques similar 
to the ones introduced in this paper and is omitted). 

From [3], we showed that a station is in need of rebalanc- 
ing if — Xi + Y^j^Li ^jPji 7^ 0- This can be easily understood 
by noting that A^ is the rate at which vehicles leave station 
i, while Y^j^Li ^jPji i s me rate at which vehicles arrive at 
station i. In what follows we assume that 

and thus each station is in need of rebalancing. We comment 



further on this assumption in Remark III. 3 

Rebalancing vehicles: In order to rebalance the number 
of vehicles Vi(t) at each station, vehicles without customers 
will be driven between stations using hired human drivers. 
The number of drivers waiting at station i is ri(t) G M>o 
and the total number of drivers in the system is R G M>o. In 
order to send a vehicle without a customer on a rebalancing 
trip from station i to station j, there must be a driver present 
at station i. We let G IR>o denote the rate at which 
we send vehicles from station i to station j when vehicles 
and drivers are available at station i. The total rate at which 



station i sends vehicles without customers is 7^ := J2j a ij> 
where an = 0. We let a denote the matrix with entries given 
by aij . These trips are shown in Figure [T] as vehicles with 
red dots in them. 

Rebalancing drivers: Finally, we must rebalance the 
drivers in the network, as they will tend to accumulate at 
some stations and become depleted at others. This is done 
as follows. If a driver would like to make a trip from station 
i to station j, it can drive a car for a customer on a trip from 
i to j, thereby acting as a taxi driver for that trip. This allows 
the driver to make the journey from station i to station j by 
"hitching a ride" on a passenger-carrying trip, but without 
negatively affecting the customer experience. We quantify 
this using two sets of variables. The variables faj G M>o 
give the rate at which drivers are sent from station i to station 
j when there are idle drivers available at station i. We let 
f3 denote the matrix with entries given by fiij and assume 

The quantities G (0, 1] give the fraction of customers 
making the trip from station i to j that would be willing 
to use the taxi mode of service on their trip. The remaining 
fraction of customers 1 — fij would prefer to drive themselves 
on their trip. Thus, imposes a constraint on the largest 
value of . In what follows we assume that the are such 
that there are enough customer trips available to rebalance 
the drivers. In Proposition |III.4 we give a necessary and 



sufficient condition on the such that this is true. These 
trips are shown in Figure [T] as vehicles with red and yellow 
dots in them. 

The notation is summarized in Table |l] 

We are now ready to write the differential equations gov- 
erning the evolution of the number of vehicles, customers, 
and drivers at each station. In order to write the expressions 
more compactly, we introduce the following notation: 

Vi := Vi(t), Ci := a(i), n := ri(i), 

:=r,(*-T, z ). q-cjQb-Tji), r) := r)(t - T^) . 

(In other words, Vj denotes the number of vehicles that were 
present at station j, specifically Tji time units prior to the 
current time.) Then, we can write the customer dynamics at 
station % as 



< 



Ai, if Vi = 0, 

0, if Vi > and q = 0, 

Xi — Hi, if Vi > and q > 0. 



Defining the Heaviside function as 
H{x) := { J; 



if x > 0, 
otherwise, 



1) the rate at which customer-carrying vehicles depart 
station i: 

-Xi, if Vi > and q = 0, 
-Hi, if Vi > and q > 0, 

which can be written more compactly as — XiH(vi) + 
{Xi-Hi)H{ci)H{vi)- 

2) the rate at which customer-carrying vehicles arrive at 
station i: 

J^Pji (A,"!'-;! - (A, - Hj)H(<j)H(v$)); 

3) the rate at which vehicles without a customer 
(rebalancing vehicles) depart station i, given by 
-jiHiv^Hin); 

4) the rate at which vehicles without a customer (re- 
balancing vehicles) arrive at station i, given by 

Thus, the vehicle dynamics can be written as 

Vi = -XiH(vi) + (Xi - Hi)H(ci)H(vi) 

+ 5>, z [x 3 H(v]) - (Xj - H 3 )H(c))H(vi)) 

- 7i H( Vi )H(ri) ■ ^n„//ir;.)//ir;.). 

Finally, the dynamics for the drivers contains four com- 
ponents. The first two components are identical to those of 
the rebalancing vehicles, given by 3) and 4) above. (This 
is due to the fact that each rebalancing vehicle contains a 
driver). The third component is the rate at which rebalanc- 
ing drivers depart station i (by driving customer carrying 
vehicles): —J2j^i^ijH(vi)H(ri). The fourth term is the 
rate at which rebalancing drivers arrive at station i with 
a customer: J2j^i PjiH( v j)H( r j)- Since drivers rebalance 
by driving vehicles on customer trips, we have from the 
customer dynamics q that 



N < 



fij^iPij if Q — 

fijfaPij if Q > 



the customer dynamics can be written as 

c % = Ai(l - H(vi)) + (Xi - Hi)H(ci)H(vi). 



However, we will consider fixed values of ftj, and since 
Hi > Xi, we simply need to enforce the more stringent 
constraint fyj < fijXiPij. 

Therefore, the fi dynamics can be written as 

h = -jiH^H^i) + Y,ajiH(v^H(r^) 



The rate of change of vehicles at station i can be written 
as the sum of four components: 



Putting everything together, we can write a set of nonlin- 
ear, time-delay differential equations describing the evolution 



of customers and vehicles in the system as 

c t =Ai(l - H(vi)) + (Xi - ^)H( Ci )H( Vl ), 
ii = - \iH{ Vi ) + (Xi - f H )H(c i )H(v i )+ 

n = - nHivjHin) + J2 o „//(,•;!//(,•;! 

where t > 0; the initial conditions satisfy q(t) = 
0, ^(r) = 0, r^r) = for r G [-max^- T ij? 0), c*(0) G 
^>o, ^i(O) G R>o with ^(0) > for at least one i G A/*, 
ri(0) G R>o with ^(0) > for at least one z G A/", 
and ^i(0) = V and Xli r i(0) = R- The optimization 
variables a and f3 are constrained as follows: 

<Pij < fij^iPij 
<otij. 

The problem we wish to solve is as follows: find an opti- 
mal vehicle rebalancing assignment a and driver rebalancing 
assignment f3 that simultaneously 

1) minimizes the number of rebalancing vehicles traveling 
in the network, 

2) minimizes the number of drivers needed, and 

3) ensures that the number of waiting customers remains 
bounded. 

Note that this is a multi-objective optimization, and thus 
it is not clear that one can both minimize the number of 
rebalancing vehicles in the network and the number of drivers 
needed. However, it will turn out that these two objectives 
are aligned, and one can find an assignment (a, /3) that 
minimizes both objectives. 

III. Well-posedness, Equilibria, and Stability of 
Fluid Model 

In this section we first discuss the well-posedness of 
model ([!} by showing two important properties, namely 
existence of solutions and invariance of the number of vehi- 
cles and rebalancing drivers along system trajectories. Then, 
we characterize the equilibria, we determine the minimum 
number of vehicles and drivers to ensure their existence, 
and we give a necessary and sufficient condition on the 
"user's preference" such that there are enough customer 
trips available to rebalance the drivers. Finally, we show 
that rebalancing vehicles and drivers give rise to equilibria 
that are locally (i.e., within a neighborhood of the nominal 
conditions) stable. 

A. Well-posedness 

The fluid model ([T]) is nonlinear, time-delayed, and the 
right-hand side is discontinuous. Due to the discontinuity, we 
need to analyze the model within the framework of Filippov 



solutions (see, e.g., [9]). The following proposition verifies 
that the fluid model is well-posed. 

Proposition III.l (Well-posedness of fluid model). For the 
fluid model ([]]), the following hold: 

1) For every initial condition, there exist continuous func- 
tions Ci(t) : R -» R>o, Vi(t) : R ->> R> , and 
ri(t) : R — >• R>o i G M, satisfying the differential 
equations ([T]) in the Filippov sense. 

2) The total number of vehicles and rebalancing drivers 
is invariant for t > and is equal, respectively, to 
^ = E^ Vi(0) andR = Y. l r*(0). 

Proof To prove the first claim, it can be checked that all 
assumptions of Theorem II- 1 in [10] for the existence of 
Filippov solutions to time-delay differential equations with 
discontinuous right-hand side are satisfied, and the claim 
follows. 

As for the second claim, the proof of the invariance of 
the number of vehicles is virtually identical to the one 
of Proposition 3.1 in [3] and is omitted in the interest 
of brevity. We prove next the invariance of the number 
of rebalancing drivers. Let r^-(t), where t > 0, be the 
number of rebalancing drivers in-transit from station i to 
station j (i.e., the rebalancing drivers for which the last 
station visited is i and the next station they will visit is 
j). Clearly, ru(t) = 0. Now, the total number R(t) of 
rebalancing drivers in the system at time t > is given, 
by definition, by R(t) = £ILi n(t) + ,. r l3 (t). The 
number of in-transit rebalancing drivers at time t is given 
by the integral over the last time units (i.e., the time 
to get from station i to station j) of the rebalancing driver 
departure rate from station i to station j. Such departure 
rate is the sum of the departure rate of rebalancing vehicles 
(since each rebalancing vehicle contains a rebalancing driver) 
and of the departure rate of rebalancing drivers that drive 
customer-carrying vehicles; hence, one can express r^ (£) as 

n 3 (t)= f aijHiviir^Hinir)) + 

J t—Tij s v ' 

rate of drivers on rebalancing vehicles 

A i ffK(T))ff(r i (r)) dr. (2) 

S ' S 

rate of drivers on customer-carrying vehicles 

By applying the Leibniz integral rule, one can write 

hjit) = K' + - H{v{)H{r{)). 

Therefore, one immediately obtains, for t > 0, 

n n n 

£(*) = EM*) + X)X>y(t) 

2 = 1 1=1 j = l 

n n 
i=l j=l 

n n n n 

i=l j = l i=l j = l 

= 0. 



This proves the claim. □ 
B. Equilibria 

The following result characterizes the equilibria of 
model ([I]). Recall that no station is exactly balanced, and 
thus -Xi + X jPji ^ °> for a11 « G 

Theorem III.2 (Existence of equilibria). Let Ax B be the 

set of assignments (a,/?) that verify the equations 

y^^ij - otji) Di, (3) 

^(Ai-fti) = -A, (4) 

/<9r £ac/z i e Af, where Di := — A^ + X^yi ^jPji- Moreover, 
let 

Vol '= ^2 Ti i (PiJ^i + a ^ 

^a,/3 := ^2 Ti i ( ai i + ^o)- 

If (a, ft) £ Ax B, then no equilibrium exists. If (a, j3) G 
Ax B, there are two cases: 

1) IfV>V a and R > R a ,p, then the set of equilibria is 

ci = 0, ^ > 0, n>0 V i G AT, 

where ^2 i v i = V - V a and = R- R a,/3- 

2) IfV<V a orR< R a ,(3, then no equilibrium exists. 

Proof. To prove the theorem, we set q = 0, V{ — 0, and 
fi = for all i e Af. From the q = equations we obtain 

A, = XiH{vi) - (A, - Hi)H{vi)H(ci). (5) 

Since \ < fa, the above equations have a solution only if 

Ci = and ^ > VieAf. 

Setting ^ = 0, combined with §5§ and the fact that in 
equilibrium q = and ^ is a positive constant, we obtain 

^2 (aijHfa) - ajiH(rj)) D u (6) 

where Di := — Xi + X^yi ^jPj*- Finally, setting = 0, 
combined with the fact that Vi > in equilibrium, we obtain 

(<XijH(n) -OLjiH^j)) 

_ (7) 

= ~ X/ i a i3 H ( r i) ~ a 3i H ( r j)) = ~ D i- 

Now, consider any station i, and note that by assumption 
we have Di ^ 0. If Di > then from ([6]) we see that Ti > 
in equilibrium. Alternatively, if Z3>^ < 0, then from ([7]) we 
see that n > 0. Therefore, in equilibrium n > 0. 

We have shown that all equilibria are of the form Ci = 0, 
^ > 0, and > 0, for each i G A/". A necessary condition 



for the existence of equilibria is that the rebalancing assign- 
ments a and (j can be chosen such that they lie in the set 
A x B of assignments that verify 

^2(0Lij - Ctji) = Di, 

I>« - &*) = ~ Di > 

for each i e Af. If (a, /3) ^ A x B, then no equilibrium 
exists and the first claim is proven. 

Assume now that (a,/?) G i x 8 and assume that 
V > V a and R > R a ,p- We need to show that c$ = 0, 
> 0, and ri > for all z G A/" are indeed valid 
equilibria. The necessary conditions in equations ^ and ^ 
are clearly satisfied and thus we simply need to verify that 
the number of vehicles and drivers are sufficient to support 
the equilibrium configuration. But, we showed in [3] that 
V a is exactly the equilibrium number of vehicles in transit. 
Similarly, from equation §2§ we can verify that R a ,p is the 
equilibrium number of drivers in transit. This, together with 
the invariance result in Theorem shows the second 

claim. 

Finally, we can show that if (a, ($) G A x B but V <V a 
or R < R a ,p, then no equilibrium exists, by arguing that in 
this case there is not a sufficient number of vehicles and/or 
drivers to support the equilibrium. □ 

Remark III.3 (Balanced stations case). We have assumed 
that Di = — Xi + ^Zj^i ^jPji 7^ for each station i. This 
assumption removes the pathological case that a station is 
perfectly balanced and does not need any rebalancing effort. 
In the case that Di = for a station, then Ti = becomes a 
valid equilibrium. Due to space constraints we have omitted 
a full treatment of the Di = case in this presentation. • 

One question remains; does there always exist an assign- 
ment (a, j3) G A x B that satisfies the constraints ol^ > 0, 
and < /3ij < fijXiPij for each z, j G A/"? We call such 
an assignment feasible. It is straightforward to verify that a 
feasible assignment for a always exists, since the variables 
are constrained only to be non-negative [3]. The /? variables, 
however, are bounded from above (that is, they have finite 
capacities), and thus it is not clear whether there exists a 
feasible (3 assignment. The following result gives a standard 
condition for the existence of a feasible assignment (see, for 
example [11, p. 220] and a consequence of this condition. 

Proposition III.4 (Existence of a feasible assignment). A 

feasible assignment (a, /3) exists if and only if, 

- ^ Dj < fij\Pij far every set S C Af, (8) 

ies ies,j<£S 

where Di = —Xi + ^j^ XjPji. As a consequence, if fij = 1 
for all Af, then a feasible assignment always exists. 

Proof. The condition ^ is a standard condition for the 
existence of a feasible solution in a minimum cost flow 
problem [11, p. 220]. 



Now we show that if = 1 for all z, j eJV, then ([5]) is 
satisfied. Take any subset S C A/" and let us show that 

A *^J Bo- 
ies' ieS,j(£S 

From the definition of .D^, the left-hand side of the above 
expression can be written as 

ies ieSjeJV iesjgs 

= ~ ^ A iPij + A -^ 

= - XiPi i + A ^ 

^2 Xipi o - °- 



of model ([]} starting at 

Ci(r) = for r G [-maxTy, 0) 

Vi(r) = z^ for r G [— maxTy, 0) 

r i( r ) = Zli f° r T ^ [— maxTy, 0) 

(c(0),v(0),r(0)) such that (1) < a(0) < ^(0) Vz, 

(2)o<r,(o) Vi,(3) 5>;(o) = v-v a , 

and (4) 5^r<(0) = R-R^ 

(11) 

has a limit which belongs to the equilibrium set. The claim 
of the theorem will then be an easy consequence of this 
statement. 

We start by observing the following fact. Assume that 
Vi(r) > and r*(r) > for all r G [- max^- T^-, t], 
then at time £ the differential equations read ti{t) = (Xi 



This proves the feasibility when = for all z, j G A/\ □ iii)H(ci(t)), for all z G A/"; recalling that, by Theorem III.2| 



it must hold 
can write 



~ x i + Ej#i A ift* - 7t + E^i «ii = 0' one 



C. Stability of Equilibria 

In this section we investigate the (local) stability of the 
equilibria of our model. We consider the following notion 
of local stability. Let (a, j3) e A x B and assume V > V a 
and > i? a ,/3 (this is a necessary and sufficient condition 



to have equilibria, see Theorem III.2). We say that the (non- 
empty) set of equilibria 



Vi(t) = -Xi + (A; - jJLi)H(Ci) - 

(X 3 - N )H(c))) 



A, 



(A, 



> ( A i - Vi)H(ci), 



for all z G A/". 



p3n 



0, ^ > 0,ri > for all 



£a,/3 •= {(c,v,r) G : 

z G A/*, and ^ = V - V a and = R- R a ,p} 

i i 

(9) 

is locally asymptotically stable if for any equilibrium 

(c, v, r) G foi 5 /3 there exists a neighborhood ^(c, v, r) := 
{(c, v, r) G R 3n | Ci > 0, vi > 0, r* > for all 'z G A/", ||(c - 
c, v — v, r — r)|| < 5, and = V — V a and = 

R — Ra,p} such that every evolution of model ([]} starting at 

c i( r ) = Q.i f° r r ^ [ — maxT^-, 0) 



Also, since by Theorem 



III.2 



it must hold - Ej^O^ 



UTone can write 



hit) 



3± 



^(t) = 2^ fc> r ^ £ [— maxT^-, 0) 



for r G [— maxT^-, 0) 



(10) 



(c(0),v(0),r(0))GB^(c,v,r)) 



has a limit which belongs to the equilibrium set. In other 
words, ( lim^ +oc c(t), lim £ ^ +oc v(t), lim t _» +00 r(t)) G 
fa,/?- The next theorem characterizes stability. 

Theorem III.5 (Stability of equilibria). Let (a, ft) G Ax B 

be a feasible assignment, and assume V > V a and R > 
Raft! then, the set of equilibria £ a ^ is locally asymptotically 
stable. 

Proof Consider an equilibrium (c, v, r) G £ a ,p (note that 
c = by Theorem |III.2| ). We now prove that every evolution 



Since Vi(r) > for all r G [— max^ T^-, 0], and since 
Vi(0) > Ci(0) for all z G AT, we conclude that no Vi(t) and 
Vi (t) can reach the value before the corresponding number 
of customers ci(t) has reached the value 0. However, once 
Ci(t) reaches the value (after a time interval Ci(0)/(fii — 
Xi)), the time derivative Vi(t) is larger than or equal to zero. 
This implies that when the initial conditions satisfy ( [TT] ), then 
Vi(t) > and n(t) > for all t > 0. 

Since Vi(t) > and Ti{t) > for all t > 0, and since 
this implies that di(t) = (Xi — Hi)H(ci(t)) for all z G M 
and t > 0, we conclude that all Ci(t) will be equal to zero 
for all t > T' := max^ Ci(0)/(fii — Xi). Then, for t > 
T' + max^j Tij =: T" the differential equations become: 
Ciit) =0, r>i(t) = 0, ri(t) =0. 

Collecting the results obtained so far, we have that 
lim £ ^ +oc Ci (t) = for all z G A/". Moreover, since ^ (t) = 
and fi(t) = for all t > T" ', the limits lim t ^+oo Vi(t) 
and lim t ^ +00 n(i) exist. Finally, one has Vi(t) = ^(0) + 
/o Vi(r) dr > Vi(0) + /* q(t) = ^(0) + a(t) - q(0). 
Since Vi(0) > q(0), we conclude that lim t ^+oo Vi(t) > 0. 
Also, fi(t) = for all t, hence lim^ +00 r^(t) > 0. Thus 
any solution with initial conditions ( [TT] ) has a limit which 



belongs to £ a ^ (the properties lim^ +00 ^Vi(t) = V — V a 
and linit^+oo ^ r *W = R ~ Ra,p are guaranteed by the 



invariance property in Proposition |IIL1| and the assumptions 
5>(0) = V - V a and ^ ri (0)=R- R^). 

Let ipi := mm(?j,, Vi sin J), and let ^ m i n := min^ 
Then, from simple a geometric argument and from the 
definitions of ipi and ^min, it follows that if one chooses 
8 — ^min, then any solution of model ([]]) with initial 
conditions satisfying fTO] ) has a limit which belongs to the 
equilibrium set. This concludes the proof. □ 

IV. Optimal Rebalancing 

Our objective is to find a rebalancing assignment (a, /3) 
that simultaneously minimizes the number of rebalancing 
vehicles traveling in the network and the number of rebalanc- 
ing drivers needed, while ensuring the existence of (locally) 
stable equilibria for model (p}. From the previous section, 
we already know that the set of assignments ensuring the 
existence of stable equilibria is A x B (provided that the 
total number of vehicles V and drivers R is large enough). 

The time-average number of rebalancing vehicles traveling 
in the network is simply given by ^j^j^j- Note that 
in minimizing this quantity we are also minimizing the 
lower bound on the necessary number of vehicles V a . The 
time-average number of drivers in the network is given by 
. Tij (oiij+ftij). Note that in minimizing this quantity we 
are minimizing the lower bound on the necessary number of 
drivers R a ,p> 

Combining the two objectives with the existence of stable 
equilibria constraints in ^ and ([4}), we obtain the following 
optimization: 

minimize ^ T^o^ and ^ (a i<7 - + /%) 



subject to y^(a^j — aji) = Di 



3& 



- Pa) = -Di 

< aij 

< Pij < fij\Pij 



VieAf 



where Di = \ + J2j^i ^jPji> anc * the optimization variables 
are and where z,j G M. The constraints ensure that 
the optimization is over the set AxB. Note, however, that this 
optimization can be decoupled into an optimization over a 
and an optimization over f3. Both optimizations are minimum 
cost flow problems [11]. The a optimization is identical to 
that presented in [3]: 



minimize 



subject to — OLji) = Di V i G J\f 

a i:j > V i, j G Af. 



The f) optimization then looks as follows: 



minimize 



subject to 



ho 

E 



-Di 



< Pij < fij\Pij 



The a optimization is an uncapacitated minimum cost flow 



problem and thus is always feasible. In Proposition III.4 
we give conditions on the fractions in order for the f3 
optimization to be feasible. 

The rebalancing policy is then given by solving the two 
minimum cost flow problems to obtain solutions a*, and 
/3*j. We then send empty rebalancing vehicles (along with 
drivers) from station i to station j at a rate of a\- (when 
vehicles and drivers are available at station i). In addition, 
we send drivers on customer-carrying vehicles from i to j at 
a rate of /?*■ (when customers and vehicles are available at 
station i). 

V. Simulations 

In this section we study the relation between the minimum 
number of drivers needed for stability R a ,p and the minimum 
number of vehicles needed V a from Theorem |III.2[ To 
evaluate these quantities, we need to generate sample data 
consisting of arrival rates at each station i, customer 
destination probabilities pij, travel times between stations 
T^, and the fraction of customers traveling from i to j 
that are willing to be driven by a driver. We generate this 
data as follows: We uniformly randomly place n stations 
in a 100 x 100 environment, and calculate the travel times 
Tij as the Euclidean distance between stations. We uniformly 
randomly generate the arrival rates on the interval [0, 0.05] 
arrivals per time unit. Similarly we uniformly randomly 
generate the destination probabilities p^j such that they are 
nonnegative and ^jPij = 1 for each station i. Finally, we 
assume that = 1 for each pair of stations in order to 
avoid issues with feasibility. 

To solve the optimizations in Section |IV| for the optimal 
assignment (a*,/?*) G A x S, we use the freely available 
SeDuMi (Self-Dual-Minimization) toolbox. 

Figure [2] shows results for numbers of stations ranging 
from 10 up to 200. For each number of stations we generate 
20 random problem instances of the form described above. 
The thick line in each plot shows the mean over the 20 
trials while the thin dashed lines show the maximum and 
minimum values. The left figure shows how V a * and R a *,p* 
vary with the number of stations. The middle figure shows 
the ratio R a *,p*/V a * as a function of the number of stations. 
We can see that we need between 1/3 and 1/4 as many 
drivers as we do vehicles. The right figure shows the ratio 
between the minimum number of rebalancing vehicles in 
transit and the number of drivers. This gives a measure of 
the fraction of drivers that are driving rebalancing vehicles 
(versus rebalancing themselves). It is interesting to note that 
this ratio is quite low, reaching approximation 1/5 for 200 
stations. 
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Fig. 2. Left figure: The minimum number of vehicles and drivers. Middle figure: The ratio between the minimum number of drivers and number of 
vehicles. Right figure: The fraction of drivers that are performing vehicle rebalancing trips. For each fixed number of stations, 20 trials were performed. 
Thick lines show the mean of the 20 trials while thin dashed lines show the maximum and minimum over the trials. 



One way to increase the fraction of drivers performing 
vehicle rebalancing is to allow multiple drivers to take 
a trip with a customer. This allows drivers to take more 
efficient routes back to stations that are in need of drivers. 
In our model it corresponds to setting > 1. This is 
explored in Figure [3] where we range from 1 to 4 for 
20 problem instances on 100 stations. We can see that as we 
increase from 1 to 4, the number of drivers decreases 
from approximately 80 to 50, and the fraction of drivers 
performing vehicle rebalancing increases from under 1 /4 to 
nearly 2/5. 

VI. Conclusions 

In this paper we studied the problem of rebalancing the 
rebalancers in a mobility-on-demand system, which blends 
customer-driven vehicles with a taxi service. For a fluid 
model of the system, we showed that the optimal rebalancing 
policy can be found as the solution of two linear programs. 



Also, we showed that in Euclidean network topologies one 
would need between 1/3 and 1/4 as many drivers as vehicles, 
and that this fraction decreases to about 1/5 if one allows 
up to 3-4 drivers to take a trip with a customer. These 
results could have an immediate impact on existing one-way 
car-sharing systems such as Car2Go. For future work we 
plan to analyze a stochastic queueing model and study the 
time- varying case whereby the system's parameters change 
periodically (thus modeling the day/night variations). Also, 
we plan to develop real-time rebalancing policies that do 
not require any a priori information, and to enrich our 
model by including uncertainty in the travel times, time 
windows for the customers, and capacity constraints for the 
roads. Finally, we are interested in using dynamic pricing to 
provide incentives for customers to perform rebalancing trips 
themselves. 




Fig. 3. Increasing the number of drivers per customer trip for 100 station problems. Left figure: The minimum number of vehicles and drivers. Right 
figure: The fraction of drivers that are performing vehicle rebalancing trips. 
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